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Particle production and its effects on tlie inflaton field are investigated during the 
preheating phase of chaotic inflation using a model consisting of a massive scalar 
inflaton field coupled to massless quantum scalar fields. The effects of spacetime 
curvature and interactions between the quantum fields are ignored. A large expan- 
sion is used to obtain a coupled set of equations including a backreaction equation 
for the classical inflaton field. Previous studies of preheating using these equations 
have been done. Here the first numerical solutions to the full set of equations are 
obtained for various values of the coupling constant and the initial amplitude of the 
inflaton field. States are chosen so that initially the backreaction effects on the in- 
flaton field are small and the mode equations for the quantum fields take the form of 
Mathieu equations. Potential problems relating to the parametric amplification of 
certain modes of the quantum fields are identified and resolved. A detailed study of 
the damping of the inflaton field is undertaken. Some predictions of previous studies 
are verified and some new results are obtained. 

PACS numbers: 04.62. +v,98. 80. Cq 

I. INTRODUCTION 

In all inflationary models there must be a period after inflation in which a substantial 
amount of particle production occurs in order to repopulate the universe with matter and 
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radiation. At the end of this process it is crucial that the universe reheats to a temperature 
that is not too large and not too small If the temperature is too large, there are problems 
with the creation of monopoles and domain walls that have not been observed, and if the 

n 

temperature is too low there can be problems with baryogenesis |2[. Thus, to assess the 
viability of a particular model of inflation, it is necessary to obtain an accurate picture of 
the amount of particle production that occurs and the rate that it occurs at. This phase of 



particle production after inflation is usually called reheating 



Both the mechanism for inflation and the way in which reheating occurs are different for 
different inflationary models. In both new inflation [sl and chaotic inflation [gI, inflation 
is due to a classical scalar field called the inflaton field. When inflation is over, this field 
oscillates in time as it effectively "rolls back and forth" in its potential 0, Q, Q]- The 
coupling of this classical field to various quantum fields and in some cases to its own quantum 
fluctuations, results in the production of particles which are thermalized by their mutual 
interactions. The backreaction of the quantum fields on the classical inflaton field causes its 
oscillations to damp j|. 

n 

The original particle production calculations for new inflation and chaotic inflation [9|] 
were done using perturbation theory. The damping of the inflaton field was taken into 
account through the use of a dissipative term that was added to its wave equation. Later, 
the mode equations for the quantum fields coupled to the inflaton field were solved and it 
was found that the time dependent part of the modes can, in some cases, undergo parametric 
amplification due to the periodic behavior of the inflaton field jiol . [ill . 12|. Since the rate 
of particle production is extremely rapid in this case and the particles do not have time 
to thermalize while it is occurring, this initial phase of the particle production process was 
dubbed "preheating" jlo| . 

The process of preheating is interesting because particle production occurs very rapidly 
with the result that quantum effects are large and the backreaction on the inflaton field is 
important. This pushes the semiclassical approximation, which is what is typically used in 
these calculations, to the limit. In fact, to take such large effects into account it is necessary 
to go beyond the ordinary loop expansion. The Hartree approximation and the large 
expansion are two ways in which such effects have been treated [l^ . 



The original calculations of preheating [1 



ll| involved solving the mode equations for 



the quantum fields in the presence of a background inflaton field which obeyed a classical 
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scalar field wave equation. Later estimates of the effects o 
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" the quantized fields on the in 
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20|. In Refs. 
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fiaton field were made using a variety of techniques 
the first fully non-linear calculations of classical infiaton decay were carried out by means 
of statistical field th eory lattice numerical simulations. The numerical simulations in 



and later ones in 23|, 
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ticeasy 



28( 1 were implemented by means of the lattice code Lat- 



291] . These simulations are important for the study of the development of equilibrium 



231 ]. Recently, the same type 



after preheating, when scattering effects cannot be neglected 
of calculation has been used to investigate the possibility that non-equilibrium dynamics, 
such as that of parametric resonance during preheating, can produce non-Gaussian density 
perturbations [30]. 

The full set of coupled equations relating the infiaton field to the quantum fields were 
solved numerically for both the Hartree approximation with a single massive scalar field with 
a quartic self-coupling and the 0{N) model using a large N expansion truncated at leading 
order. The 0{N) model consists of N identical massive scalar fields each with a quartic 
self-coupling and a quadratic coupling to every other field. As a result there are N infiaton 
fields in this model. The first set of calculations were done in a Minkowski spacetime back- 



ground 
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3l| and then followed by radiation and matter dominated Friedmann-Robertson- 



Walker backgrounds [32!]. As the fields have an energy- momentum tensor associated with 
them, the expansion of the universe is also affected. This was taken into account in Ref. 
A numerical calculation in the context of the SU(2) Higgs model was done in Ref. 3j]. A 
calculation within the context of the QiN) model that takes scattering of the produced par- 
ticles into account was done in Ref. 35|]. Numerical calculations in models with two scalar 
fields, with backreaction effects taken into account using the Hartree approximation, were 
done in Ref. [36 ]. 

There are some potential problems with using the 0{N) model as a model for chaotic 
inflation. One is that it results in infiaton fields rather than the single infiaton field 
that is usually postulated. Another is that the fields have a quartic self-coupling as well 
as couplings to each other. There is evidence from the WMAP III data [37] that if chaotic 
inflation occurred, the infiaton field is a massive scalar field with no quartic self-coupling. 
Finally, in the 0{N) model or in a model in which the primary quantum fluctuations come 
from the infiaton field, the quantum fields are massive with the same mass as the infiaton 
field. However, in a realistic situation there would also be massless or effectively massless 
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fields coupled to the infiaton field. 

A possibly more realistic model for chaotic inflation consists of a single classical infiaton 
field coupled to various quantum fields, including its own quantum fluctuations. This is 

nnn 

the type of model considered by Kofman, Linde, and Starobinsky p^, [UJ, 115[ . They used 
solutions of the mode equations of the quantum fields to estimate both the amount of particle 
production that would occur in various cases and the amount of damping of the infiaton 
field that occurs due to the backreaction of the quantum fields on it. 

In this paper we consider a model similar to one used by Kofman, Linde, and Starobinsky, 
henceforth referred to as KLS, in their analysis of preheating in chaotic inflation [l^. The 
infiaton field is a classical massive scalar field with minimal coupling to the scalar curvature. 
It has no quartic or higher order self-coupling. But it is coupled to identical massless 
quantized scalar fields with arbitrary curvature coupling C,- As discussed in Section II, in 
the large A^ limit, the system effectively reduces to one consisting of a classical infiaton field 
coupled to one quantized massless scalar field with arbitrary curvature coupling. The effects 
of quantum fluctuations of the infiaton field come in at next to leading order so we do not 
take them into account. 

We investigate some of the details of the preheating process by numerically integrating 
both the mode equations for the quantized field and the backreaction equation for the 
infiaton field. These are the first solutions to the full set of coupled equations relating the 
infiaton field to the quantum field that have been obtained for this model. To simplify 
the calculations we ignore interactions between the created particles. Calculations that 



have taken such interactions into account 
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351 ] indicate that they can be 



ignored during the first stages of preheating but eventually become important. Another 
simplification is that, as a first step, we work in a Minkowski spacetime background. It 
has previously been shown that the expansion of the universe can significantly affect the 
evolution of the infiaton field during preheating if this field is of Planck scale or larger at the 
onset of inflation [l^ . This in turn can substantially affect the details of the particle creation 
process. Thus, some of our results will not be relevant for most models of chaotic inflation. 
However, for the rapid damping phase which occurs in many models, the amplitude of the 
inflaton fleld changes rapidly on time scales which are small compared to the Hubble time. 



Thus, it should be possible to ignore the expansion of the universe during this phase 



15|- 



If, for the model we are considering, one neglects the backreaction effects of the quantized 
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fields on the inflaton field, then the inflaton field undergoes simple harmonic motion with a 
frequency equal to its mass. As shown in Section III, this means that the time dependent part 
of the modes of the quantized field obey a Mathieu equation. Thus, parametric amplification 
will occur for modes in certain energy bands. There is an infinite number of these bands 
and some occur for arbitrarily large values of the energy. This could in principle result 
in divergences of quantities such as (V^^), with ip the quantum field. For example, if the 
contribution to the mode integral in {ip"^) [see Eq. (12. lip ] for each band of modes undergoing 
parametric amplification was the same at a given time, then {ip"^) would diverge. Another 
issue that must be addressed is that even if {ip"^) is finite, one must be certain that all 
significant contributions to it are accounted for when making numerical computations; in 
other words, one must be certain that no important bands are missed. These issues are 
addressed in Section III. 

When the backreaction of the quantized fields is included in the wave equation for the 
inflaton field, then, as expected, the inflaton field's amplitude damps as the amplitudes of 
some of the modes of the quantized field grow. KLS predicted that if all the instability 
bands are narrow then the damping is relatively slow while if one or more of the instability 
bands are wide then there is a period of rapid damping. We find this to be correct, and 
for a Minkowski spacetime background, we find a more precise criterion which determines 
whether a phase of rapid damping will occur. In cases where there is rapid damping, we find 
a second criterion that must be satisfied before it takes place. In a Minkowski spacetime 
background, this latter criterion explains the differences that are observed in the time that 
it takes for the period of rapid damping to begin. 

In their study of the rapid damping phase, KLS gave a criterion for when the damping 
should cease and used it to predict how much damping should occur. We find something 
similar except that the amount of damping that is observed to occur is larger than they 
predict. We also find that the rapid damping actually occurs in two phases separated by a 
short time. An explanation for this is provided in Section IIVDI 

If most of the damping occurs gradually then it is observed that the frequency of the 
oscillations of the inflaton field changes slowly as its amplitude is damped. There is also a 
significant transfer of energy away from the inflaton field as would be expected. However, 
if most of the damping occurs rapidly then during the rapid damping phase the frequency 
of the oscillations of the inflaton field increases significantly. As a result there is less energy 
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permanently transferred away from the inflaton field than might be otherwise expectec . 



22|. 



This result was seen previously in the classical lattice simulation of Prokopec and Roos 

No significant further damping of the amplitude of the inflaton field was observed to occur 
after the rapid damping phase in all cases in which such a phase occurs. Instead both the 
amplitude and frequency of the inflaton field were observed to undergo periodic modulations 
and a significant amount of energy was continually transferred away from and then back to 
the inflaton field. However the actual evolution of the inflaton field after the rapid damping 
phase is likely to be very different because interactions that we are neglecting are expected 



to be important during this period [2l|, l22|, l23|, |26|, |27|, 1351 ]. 



In Section II the details of our model are given and the coupled equations governing the 
inflaton field and the quantum fields are derived. In Section IIIII a detailed study is made 
of parametric amplification and the bands that undergo parametric amplification in order 
to address the issues discussed above relating to the finiteness of certain quantities and 
the potential accuracy of numerical computations of these quantities. In Section IV some 
of our numerical solutions to the coupled equations governing the modes of the quantum 
field and the behavior of the inflaton field are presented and discussed. Our results are 
summarized in Section V. In the Appendix we provide the details of the renormalization and 
covariant conservation of the energy-momentum tensor for the system under consideration. 
Throughout this paper we use units such that h = c = G = 1. The metric signature is 



II. THE MODEL 



We consider a single inflaton field $ and identical scalar fields \E'j, that represent the 
quantum matter fields present during the inflationary phase. We assume that the inflaton 
field is massive (with mass m) and minimally coupled to the scalar curvature, R, and that 
the quantum scalar fields are massless with arbitrary coupling C, to the scalar curvature. 
We also assume that the interaction between the classical inflaton and the quantum matter 
fields is given by (7^$^\E'|/2. We study the dynamics of $ and in a Minkowski spacetime 
background and denote the metric by rjn^- 
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The action for the system is given by 

1 



N 



(2.1) 



One can divide the fields into classical and quantum parts by writing $ = 0c + and = ijj 
with ($) = (f)c and = 0. Then if the infiaton field is rescaled so that 0c — > \/N(f)r it is 
possible to carry out a large N expansion^ using the closed time path formalism 39|. The 
result to leading order is that the system is equivalent to a two field system with the classical 
infiaton field (pc coupled to a single quantized field ip, henceforth referred to simply as (p and 
ip, respectively. The equations of motion for these fields are 



(-□ + m2 + /(^2^B)0 = 



(2.2a) 
(2.2b) 



where {iP'^)b is the bare (unrenormalized) expectation value of ip"^. Restricting to the case 
of a homogeneous classical infiaton field in Minkowski spacetime, the equation of motion 
simplifies to 



(Pit) + {m' + g'{r)B)(l>it) = 0. 



(2.3) 



Equation ( ]2.2b[) is separable and the quantum field tpi^) can then be expanded as 

d^k 



V^(x,t) 



ake^'^-^/,(t) + aj,e-^^V:W 



(27r)3/2 

The time dependent modes, fk{t), satisfy the following ordinary differential equation 

fk{t) + [k' + g'^'mkit) = 0, 



and the Wronskian (normalization) condition 



(2.4) 



(2.5) 



(2.6) 



^ The expansion here is similar in nature to one that was carried out for Quantum Electrodynamics with 
N fermion fields in Ref. 
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The WKB approximation for these modes is useful both for renormahzation and fixing 
the state of the field. It is obtained by writing 



fk{t) 



1 



exp 



v/2Wfc(t) 

Substitution of Eq. ([23D into Eq. gives 



2 i2 



-I I dt' Wkit') 

to 



3W, 



(2.7) 




2 



(2.8) 



One solves this equation by iteration. Upon each iteration one obtains a WKB approxima- 
tion which is higher by two orders than the previous one. Note that the order depends both 
upon the number of time derivatives and the power of g. The zeroth order approximation 



is just Wf 



(0) 



k, while the second order approximation is 

,2 J,2 



(2) 



k + 



g 



2k 



(2.9) 



and the fourth order one is 



(4) 



k + 



g 



2^2 



g 



4:M 



3_ 

4P 



+ 



2k 8A;3 

Eq. (12. 3p for the inflaton field involves the quantity {ip'^)-B- Using (12.41) one finds 

1 



(2.10) 



{r)i 



27r2 



dkk' \h{t)[ 



(2.11) 



This quantity; 



ularization 



40 



is divergent and must be regularized. We use the method of adiabatic reg- 



41 



42 



43| which, for free scalar fields in Robert son- Walker spacetimes has 



been shown to be equivalent to the covariant scheme of point splitting 4J, |45l] 



In adiabatic regularization the renormahzation counterterms are o 



3tained by using a 



WKB approximation for the modes of the quantized field |42|, |43|, |4J, |45|, |46|, |47|. One 
works with a massive field and then takes the zero mass limit at the end of the calculation. 
Expressions for the renormalized values of both {ip"^) and the expectation value of the energy- 
momentum tensor {T^J) for a similar system to the one we are using^ have been obtained 



^ In Ref. ^47,], if the scale factor, o(t), is set equal to one, the coupling of the inflaton field to the scalar 
curvature is set to zero, the coupling constant A is set equal to 2g^, the mass of the quantum field is set 
to zero, and the A0'^/3! term in the equation of motion for the (f) field is dropped, then the two systems 
are equivalent. 
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in Ref. 47|. Some details of the renormalization procedure are given in the Appendix. The 
result for our system is 



0^% 



27r2 



,2 J,2 



+ 00 



,2 J,2 



2k 4A;3 



(2.12a) 



(2.12b) 



where e is a lower limit cutoff that is placed in the integrals that are infrared divergent 
and M is an arbitrary parameter with dimensions of mass which typically appears when the 
massless limit is taken 



46 



47l | . In principle the value of M should be fixed by observations. 
For simplicity, we set it equal to the mass of the inflaton field, m, in the calculations below. 
Note that the value of (V'^)ren is actually independent of the value of the infrared cutoff e. 

It is useful for both the analytic and numerical calculations to scale the mass of the 
inflaton field, m, out of the problem. This can be done by defining new dimensionless 
variables as follows 

(2.13a) 
(2.13b) 

(2.13c) 
(2.13d) 
(2.13e) 
(2.13f) 

If one substitutes (V'^)ren for {iIj^)b in Eq. (EJl), writes Eqs. (EJ]), (ESI), (EBl), and (KWf 
in terms of the dimensionless variables, rescales {ip^)ren so that {ip'^)ren "^~^(V'^)ren5 and 
then drops the bars, the equation for the inflaton field [Eq. (12. 3p ] becomes 



t = 


mt , 


= 


m ' 




k 


k = 






m ' 


fk = 






M 


M = 






m ' 




e 


e = 






m 



0(t) + [l+^72(^2^,en]0(t) 







(2.14) 



while Eqs. (12. 5p . (12. 6p and (12.120 remain the same. These four equations in terms of the 
dimensionless variables are the ones that are solved in the following sections. 
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III. BACKGROUND FIELD APPROXIMATION 

In the background field approximation the wave equation for the inflaton field is solved 
without taking into account the backreaction effects of the quantum fields on the inflaton 
field. However, the effect of the inflaton field on the mode equations for the quantum fields 
is taken into account. For our model the wave equation for the inflaton field (12.141) then is a 
simple harmonic oscillator equation, and the mode equation Eq. (12.51) is a Mathieu equation. 

It is well known that for the Mathieu equation there are regions in parameter space 
for which there are solutions which grow exponentially due to a process called parametric 
amplification. For the mode equation (12. 5p the last term is proportional to g'^ct)^., with 0o the 
amplitude of the oscillations of the inflaton field. Thus, for a given (^^^q, it is the parameter 
k which determines which modes undergo parametric amplification. There are bands of 
values of k for which this occurs and we shall call them instability bands. The instability 
band which contains the smallest values of k will be called the first band. 

For a given value of g'^cf^ there is an infinite number of instability bands some of which 
contain modes with arbitrarily large values of k. In principle, this could lead to divergences 
in quantities such as (V'^)ren and the expectation value of the energy-momentum tensor, 
(T^j,)ren- The rcasou is that, as can be seen from Eq. (I2.12al) . the renormalization countert- 
erms do not undergo any exponential growth while bands of modes with arbitrarily large 
values of k do. One might also be concerned that, even if (V'^)ren and (T^i/)ren are are finite, 
in order to compute them one must take into account modes in instability bands with ar- 
bitrarily large values of k. This is clearly not possible for purely numerical computations, 
which must necessarily include only a finite number of modes and for which there must be 
an ultraviolet cutoff. In this section we use known properties of solutions to the Mathieu 
equation to show that both {4'^)ren and {T^u)ren are finite at any finite time. We also address 
the question of how to make sure that the contributions of all of the instability bands which 
contribute significantly to (V'^)ren and (T^:y)ren are included in numerical calculations of these 
quantities. 



11 



A. The Mathieu Equation 



In the background field approximation, without loss of generality (due to time translation 
invariance), we can take the solution to Eq. f l2.14p to be 



= </)ocos(t) . 



The mode equation is then 



Using the identity 



1 



cos^(t) = -[l + cos(2t)] 



the mode equation can be put into a standard form for the Mathieu equation 



+ [a-2qcos{2t)]fk{t) = 



with 



; 2 I 2 J 2 



2 i2 

Q = -^9 00 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



(3.5a) 
(3.5b) 



Floquet's theorem applies to the Mathieu equation 49|]. It implies the existence of two 
solutions of Eq. (13. 4p of the form: 



m 

f2{t) 



(3.6) 
(3.7) 



where hAt_ 



number 



and h2{t) are periodic functions of period 27r and /x is, in general, a complex 
491]. If has a no n- vanishing real part, it is easy to see that one solution will grow 



exponentially at late times. 

Our goal is to find which ranges of values of k lead to exponentially growing modes. 
This can be accomplished by first fixing the value of q and then finding the values of the 
parameter a for which the Mathieu equation has solutions with /i = 0. Those values of a 



12 



which correspond to even solutions are traditionally labeled and those which correspond 
to odd solutions are labeled b^, with r a positive integer. The pattern of instability and 
stability regions in the (a, q) plane is symmetric about q = 0. Thus, for simplicity we shall 
consider g > in what follows, even though g < for the equations we are concerned with 
[see Eq. (I3.5bl) ]. For q > O^jt is found that /x^ > for values of a which are between and 



br for the same value of r [48||. Further one has the relation oq < &i < Oi < 62 < • • • • 

For the purpose of assessing the contribution of unstable modes with large values of k to 
Eq. (12.121) . it is necessary to know both the range of values of k in the unstable bands and 
the maximum value of /x in a band containing large values of k. The ranges of k in the large 



r limit are given by the relation 50] 



ar-br = . (3.8) 

Thus, for fixed q, the width of the unstable bands in the (a, q) plane becomes arbitrarily 
small for bands with large values of r and therefore becomes arbitrarily small for large values 
of k. This means that for unstable modes at large values of k to contribute significantly to 
the integrals in Eq. (12.121) . it would be necessary that these modes grow much faster than 
the unstable modes in bands with smaller values of k, which are wider. 



An approximation to the value of fi in unstable bands is 48] 



2 {ar + br-2a- Ar"^) ^ al - 2aX - ^r\ar + br - 2a) + 6^ + iQr^ 
/^r - 2 2 ' 

where the boundaries of the unstable band of index r are given by the pair of eigenvalues br 
and Oj, with > br- Some details of the derivation of this approximation can be found in 
Ref . [511 • Note that there are two possible solutions for /i^ for fixed values of r and a. Given 
that for large r the instability bands become narrower, and that br < a < ar, it is easy to 
show that for large r only the plus sign gives yU^ > 0. It is also easy to show that 

\2a - ar - br\ 4:r^ , (3.10) 

and thus that 

f^r - . (3.11) 

For this approximation the maximum value of fir is at the midpoint between br and ar- 
Since, for large r the width of the band, a.^ — br, gets very small [see Eq. (13.81) ]. it is clear 
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that the maximum value of /i in a band is smaller for bands with very large values of k than 
for bands with smaller values of k. 

The relative contribution of unstable modes to the integral in Eq. f l2.12p at a given time t 
depends in part on how much they have grown compared to unstable modes in lower bands. 
It is clear from the above analysis that modes in higher bands grow significantly slower than 
those in lower bands. Also for a given time t, there will be some value, k = K, for which 
unstable bands with larger values of k have not grown much at all. This, coupled with the 
fact that the widths of the unstable bands get narrower at larger values of k, means that 
at any given time t the modes in bands with k > K will not contribute significantly to the 
mode integral in Eq. (12.121) . Thus, (V'^)ren and {Tfj_u)ren are finite. 

From the above analysis it should be clear that in numerical computations it is possible 
to obtain a good approximation to {ip'^)rcn and (T^,y)rcn by putting in a cutoff at some large 
value of k. Since the location of the bands at small values of k can be determined using 
standard techniques [481], then one simply needs to use a high enough density of modes so 
that those bands with smaller values of k are adequately covered. 

The maximum value of fi occurs for the first unstable band which has the smallest values 
of k. Once the modes in this band have become sufficiently large, they make the major 
contribution to (V^^)ren and {T^u)Tcn- 



IV. BACKREACTION ON THE INFLATON FIELD 



A. Initial Conditions 



A Fortran program was written to simultaneously solve the mode equations for the ip field, 
compute the quantity {ip'^)r:en, and solve the equation for the infiaton field in terms of the 
dimensionless variables (12.131) . For this program to run, initial values must be given for (p and 
its first time derivative, as well as each mode function and its first time derivative. Since 
the infiaton field is in its rapid oscillation phase, what is important is the initial amplitude 
of its oscillations. This is most easily obtained by starting with 0(0) = 0o and 0(0) = 0. 

The initial values for the mode functions of the ip field are determined by the state of 
the field. As the quantum field is coupled to the infiaton field, there is no preferred vacuum 
state. Instead one can use an adiabatic vacuum state. Such states have been discussed 
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in detail for cosmological spacetimes [52j. They are based on the WKB approximation for 
the modes. It is necessary to choose a state that is of adiabatic order two or higher for 
the renormahzed value of {ip'^) to be finite. For the renormalized energy- momentum tensor 
(T^jy) to be finite it is necessary to have at least a fourth order adiabatic state. For such a 
state, in the large k limit, 

|.p^±_i^^ + 5i^*! + 9!«M + o(t-») (41) 

'•''=' 2* ^ i6*!> ^ se ^ ' * ' 

Although the large k behavior of an adiabatic state of a given order is constrained by 
requiring that renormalization of the energy-momentum tensor be possible, there are no 
such constraints on modes with smaller values of k. Thus, one can construct fourth order 
adiabatic states that would give virtually any value for {ip'^)ren or (T^jy)rcn at some initial 
time t. However, if the value of g'^{ip'^)ren is too large at the beginning of preheating, 
then backreaction effects will be important immediately. In fact, if g'^{ip'^)ren ^ 1) then 
backreaction effects are so strong initially that parametric amplification is unlikely to occur. 
If g"^ {ip"^) icn ~ 1, backreaction effects will be important, but it is possible that parametric 
amplification might still occur. On the other hand if 5'^(V'^)ren <^ 1 at the onset of preheating, 
then backreaction effects are not important initially and parametric amplification will occur. 
Since our goal is to study the preheating process in detail, we restrict attention to this latter 
case. 

A natural way to construct a fourth order adiabatic state is by using a fourth order 
WKB approximation to fix the initial values of the modes, fk, and their derivatives, fk, for 
all values of k. For the numerical calculations done in this paper the state used was an 
adiabatic vacuum state with 

fk{t = 0) = (4.2a) 
Mt = 0) = --^ , (4.2b) 

Y = — + (4 2c) 



^^0 = ^Jk^ + g^l. (4.2d) 

It is easy to show that this is a fourth order adiabatic state by recalling that 0(0) = 0. The 
initial values have been chosen in part to make it possible to analytically compute the initial 
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value of ('?/'^)ren5 and in part so that the Wronskian condition (12.60 is satisfied exactly. For 
the fourth order adiabatic state that we are using, at the initial time t = 



log 



M2 



60 levr^ 

Using this initial value for 5'^(V'^)ren one finds that the initial value for is 

1 



(4.3) 



0(t = 0) 



1 + 



487r2 



-1 + ^11 

167r2 



log 



M2 



(4.4) 



Thus, whether backreaction effects are important initially depends on the values of both 0o 

2 

and g. For the values of g considered here, <^ 1. From Eq. (14.41) it can be seen that 
backreaction effects for the states we are considering are likely to be important if 



g\g^<Pl) > levr^ . 



(4.5) 



As an example, for realistic models of chaotic inflation KLS 15[ chose 



m 



5m 



(4.6) 



400, while 



The values they used for g range from ~ lO^'^ to 10^^. For g = 10~^, g'^(po 
for g = 10~^, g'^(pQ = 4 X 10^. Applying the criterion (14. 5 p one flnds that our fourth order 
adiabatic states result in signiflcant backreaction initially for these values of m and (j)^ if 
^ > 8 X 10^3_ 



B. Parametric Amplification With Backreaction 

If the initial conditions are such that 5'^(V'^)ren 1, then examination of Eq. (12.141) shows 
that there is no signiflcant backreaction at early times and the solution of the inflaton fleld 
will be approximately equal to Eq. (13.11) . When, primarily through parametric ampliflcation 
of certain modes, g'^{ip'^)Ten ~ 1, then one expects backreaction effects to be important and to 
make the amplitude of the oscillations of the inflaton fleld, 0, decrease. It has been observed 
in the numerical computations that {ip'^)rcn oscillates about an average value which changes 
in time. The effective mass of the inflaton fleld can be deflned in terms of this average value: 

mis = 1 + g'm.en ■ (4.7) 
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It is clear that an increase in results in an increase of the frequency of oscillations of 
(f). This should happen slowly at first, and so one would expect the basic character of the 
mode equation fl2.5p to be a Mathieu equation, except that the effective value of q would 
change 15|. A zeroth order WKB analysis of the solutions to Eq. (12.141) shows that the 
effective frequency of oscillations can be approximated by 

1 



t2 

dt v^TTTPU , (4.8) 



h — h jt^ 

with ti < t < t2 for some time interval t2 — 1\ equal to or larger than an oscillation period of 
the inflaton field. A derivation similar to that at the beginning of Sec. Ill then shows that 

fleff ~ — + 2geff , (4.9a) 

'^efT 

where A is the amplitude of the inflaton field and we have defined gefr so that it is positive^. 
As the inflaton field evolves A will tend to slowly decrease and cUcfr will tend to slowly 
increase. Thus, the resulting changes in the instability bands should be small enough that 
there would be no problem in having enough modes to accurately determine the primary 
contributions to (V'^)ren- Our numerical work described below seems to bear this out. 

Once the backreaction is significant, it is possible for the character of the equation to 
change substantially. If varies in a periodic manner then it is still possible for parametric 
amplification to occur due to Floquet's Theorem 49|. It would be difficult in this more 
complicated situation to analyze the behavior of the solutions to the extent that solutions 
to the Mathieu equation have been analyzed. However, it seems quite likely that, as for the 
Mathieu equation, the most important instability bands will be at relatively small values of 
k and that if an instability band is extremely narrow, then it probably will not have a large 
effect on ('?/'^)ren for a very long time. It is likely that interactions which are neglected in 



our model, will be important on such long time scales |2lL 1221 . l23l . 1261 . l27l . ISSj 



y. 



^ The actual derivation results in a negative value for gcff but the symmetry between positive and negative 
values of q for the Mathieu equation allows us to simply consider positive values. 
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C. Some Numerical Results 

Since the mode equation (12. 5p depends on g^cfp' and the equation of motion of the inflaton 
field f l2.14p depends on g'^{ip'^)^en, one would expect that the effects of the values of g and 
5f^0o decouple, at least to some extent. As can be seen in Figs. [H and |2] this is the case, 
with the quantity g^(pl determining the type of evolution the inflaton field undergoes and g 
having its strongest effect on the timescale of that evolution. 

In Fig. [1] plots of solutions to the backreaction equation fl2.14p for the case in which 
g = 10~^ are shown, with g'^(pl ranging from 35 to 1. For the first two plots most of 
the damping occurs over a very short period. After that, the field does not appear to 
damp significantly. Instead, it continues to oscillate but with a much higher frequency than 
before. The envelope of its oscillations also oscillates but with a much smaller frequency. 
This behavior has been observed for all cases investigated with g'^(f)Q ~ 2. Conversely, for the 
plot on the bottom {g'^(f)l = 1), it is clear that the damping is much slower and occurs over 
a much longer period of time. This behavior is a verification of the prediction by KLS [l^ 
that a period of rapid damping would occur for g'^(f)Q ^ 1 and that only relatively slow 
damping occurs for g^cj)^ <^ 1. A detailed analysis of the evolution of the inflaton field is 
given in subsection lIVDi 

In Fig. [2] plots of solutions to the backreaction equation for the case in which g'^(f)l = 10 
and g = 10~^ and 10~^ are shown. The amount of damping does not appear to depend in 
any strong way upon the value of the coupling constant g. The amount of time it takes for 
significant damping to begin is longer for the smaller value of g. This is easily explained 
by the fact that initially g'^{ip'^)ren ^ 1? so that only negligible damping of the inflaton field 
occurs at early times. During this period the mode equations and thus the growth of (^/'^)ren 
is approximately independent of the value of g. Significant backreaction begins to occur 
when g"^ {ip"^) ren ~ 1 and this will obviously take longer for smaller values of g, since (V'^)ren 
will have to grow larger before it can occur. 

The natural explanation for what is happening in all of these cases is that particle pro- 
duction occurs and the backreaction of the produced particles on the inflaton field causes 
its amplitude to be damped. However, if this picture is correct then one would expect that 
a significant amount of energy would be transferred from the inflaton field to the modes of 
the quantum field. To check this, one can compute the time evolution of the energy den- 
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FIG. 1: The evolution of the inflaton field for g = 10 ^. From left to right the top plots are for 



35 and 10. The bottom plot is for ^^^n ~ 



sity of the inflaton field. Since the gravitational background is Minkowski spacetime, the 
energy density is conserved as is shown explicitly in the Appendix. As is also shown in the 
Appendix, the energy density can be broken into two different contributions. The first one 
is the energy density the inflaton field would have if there was no interaction (^f = 0). It is 
given by 

(4.10) 



()2 + mV 



The other contribution is the energy density of the quantum field shown in Eq. (]A6cp . 
It contains terms that would be there if there was no interaction along with terms that 
explicitly depend on the interaction. This split is useful because almost all of the energy is 
initially in the inflaton field and thus, one can clearly see how much has been transferred to 
the quantum field ip and to the interaction between the two fields as time goes on. 

In Fig. [3] the left hand plot shows the evolution of the energy density of the inflaton field 
for the case g = 10^^ and g'^4'Q = 1. Comparison with the bottom plot in Fig. [1] shows that 
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FIG. 2: The evolution of the inflaton field for 5 = 10- The plot on the left is for for g = 10 



and the one on the right is for 10 ^. 
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FIG. 3: The evolution of the energy density of the inflaton field, Pfj,. Both plots are for g = 10 ^. 
The one on the left is for g'^cpQ = 1 and the one on the right is for g'^cpQ = 10. 

as the inflaton field is damped, a large portion of its energy is permanently transferred away 
from it. 

For all cases investigated in which rapid damping of the inflaton field occurs, it was found 
that much less energy is permanently transferred away from the inflaton fleld than one might 
expect, given the amount by which its amplitude has been damped. This effect was seen 
previously in the classical lattice simulation of Prokopec and Roos 22] . An example is the 



case g = 10^^, g'^cf^ = 10 for which the evolution of the energy density of the inflaton fleld is 
shown in the plot on the right panel of Fig. [3l Comparison with the top right plot of Fig. [2] 
shows that after the rapid damping has flnished the inflaton fleld permanently loses some 
energy because the maximum of its energy density is lower than before, but more than half 
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FIG. 4: This figure shows differences in (/> and g^(V'^)ren which occur for g = 10-3 and c/^^g ^ iq 
and 5'^</>o = 10(1 + 10^^). In each plot the dashed curve corresponds to g^(l)Q = 10(1 + 10^^). 

of the original energy is transferred back and forth between the inflaton field, the quantum 
field, and the interaction between the two. However, it is important to note that in cases 
in which rapid damping occurs our results cannot be trusted after the rapid damping phase 
has finished because, as previously mentioned, there is evidence that the interactions we 



neglect become important 
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Another effect that was observed to occur in all cases in which the inflaton field undergoes 
rapid damping, is an extreme sensitivity to initial conditions. If the initial conditions are 
changed by a small amount then initially the evolution of the inflaton field and the quantity 
('?/'^)ren are nearly the same as before. However, at some point the evolution of {ip'^)ren begins 
to change significantly and that is followed by a significant change in the evolution of 0. An 
example is shown in Fig. HI This phenomenon is discussed in more detail below. 



D. Detailed Analysis 

As discussed in Section lIVBt as {i''^)ren increases due to parametric amplification the 
amplitude of the field damps and its oscillation frequency increases. Thus, it is useful to 
define an effective value for |g|, which we call qes, and whose definition is given in Eq. (14.9bl) . 
If Qes ^ > for some or all of the modes in an instability band, the oscillation frequency 
for the modes is much larger than that for the inflaton field. KLS point out that this allows 
one to use the Mathieu equation to provide a reasonable description of what happens since 
geff varies slowly compared to the oscillation time for the modes. If ges is small enough 
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FIG. 5: This figure shows the boundaries of the five lowest instabihty bands for the Mathieu 
equation using (j3.5ap and ()3.5bp . The lower bound of a given band is the dashed curve and the 
upper bound is the solid curve immediately above it. Note that there are no modes of the quantum 
field with k"^ < 0. 

then it seems likely that a description in terms of the Mathieu equation breaks down from 
a quantitative, and perhaps also from a qualitative, point of view. 

Our explanations of the details of the damping process are based on the basic scenario 
outlined by KLS in which as the damping proceeds, qes decreases and there are various 
instability bands that cross over a given set of modes. The damping stops when gefr gets 
small enough. The locations of the five lowest instability bands in the {k"^, q) plane are shown 
in Fig. [51 

There are four primary factors that go into our explanations: 

1. The lowest instability band, i.e., that with the smallest values of k, has the largest 
value of /imax and is the widest. Thus the modes near its center undergo the largest 
amount of parametric amplification. 

2. As Qcs decreases the tendency is for instability bands to move to the right along the k 
axis and for new ones to appear when k^ = u'^g^acs — 2qes) becomes positive for part 
of the band. 

3. Modes near the center of an instability band will continue to undergo significant para- 
metric amplification (compared with the rest of the band) until the backreaction on 
the infiaton field causes gefr to decrease enough that these modes are either near the 
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FIG. 6: These plots show (7^(V'^)ren for g = 10^^ and 5'^</>o = 1- The plots are identical except that 
the one on the left is plotted on a logarithmic scale and the one on the right is plotted on a linear 
scale. 

edge of the instability band (where /i is small) or no longer in it, due to the shift of 
the band along the k axis. 

4. Once the lowest instability band shifts significantly so that the original modes near 
the center stop growing rapidly, new modes will begin to undergo a significant amount 
of parametric amplification. However, it takes a while for the contributions of these 
new modes to the value of {ip'^)ren to become significant compared with that of the 
modes which were previously near the center of the first instability band, and have 
now stopped undergoing significant parametric amplification. 

In what follows we first discuss the case in which no significant amount of rapid damping 
occurs. For the cases in which rapid damping does occur the pre-rapid damping, rapid 
damping, and post-rapid damping phases are discussed separately. 

KLS predict and we observe that there is no significant amount of rapid damping for 
^ 1. However, the details are different because they consider an expanding universe. In 
the Minkowski spacetime case, as can be seen by comparing the lower plot of Fig. 1 and the 
plots in Fig. E], there is actually a small amount of rapid damping that occurs for g^cpl = 1 
once g"^ {ip"^) ren ~ 10~^. What appears to be happening is that modes near the center of 
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FIG. 7: This figure shows the integrand in Eq. (j2.12ap for the case g = 10~^, g'^cpQ = 1, for the 
times t = 150 (left) and t = 250 (right). Note that at the later time the band of modes that 
contributes significantly to the integral includes modes with larger values of k than the band at 
the earlier time. 

the first instability band cause ('?/'^)ren to grow exponentially until the band has shifted to 
the right along the k axis enough so these modes no longer grow so rapidly. Because of the 
exponential growth of (7^('?/'^)ren, the damping of the infiaton field is very rapid once it begins 
to occur in a significant way. However, as |g| = 0.25 in this case, the first instability band 
is relatively narrow and thus not much damping of the infiaton field needs to occur for this 
band to shift by the required amount. Once it has shifted, new modes begin to undergo 
significant parametric amplification and as they start to contribute to {ip'^)rcn the amplitude 
of the infiaton field continues to decrease and the position of the instability band continues 
to move slowly to the right along the k axis. One can see evidence of this shift by comparing 
plots of the integrand for {ip'^)rcn at two different times as is done in Fig. [71 

As the infiaton field continues to damp and ges continues to decrease, the value of fi at 
the center of the band should also be decreasing. Thus, the time scale for the new modes 
undergoing parametric amplification to contribute significantly to the damping increases and 
the rate of damping decreases. 

2. g''(g>2 

For (7^00 ~ 2 rapid damping has occurred in every case that has been numerically inves- 
tigated. However, the time at which the rapid damping occurs and the details relating to 
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how the rapid damping occurs and how much damping there is vary widely. To understand 
the reasons for this, it is useful to break the discussion up into the pre-rapid damping, rapid 
damping and post-rapid damping phases. 

a. Pre-Rapid Damping Phase 

Since g'^{ip'^)ren ^ 1 initially for all cases considered in this paper, there is no significant 
amount of backreaction on the inflaton field at first. Therefore, parametric amplification of 
the modes near the center of the first instability band, i.e., the one which encompasses the 
smallest values of fc, occurs and the contribution of these modes quickly dominates all other 
contributions to {ip'^)ren- ^ (As mentioned previously this is in contrast to the case of an 
expanding universe where the expansion causes a decrease in the amplitude of the inflaton 
field and thus in gofr-) As before, (V'^)ren grows exponentially until it becomes of order 10~^ 
at which point significant damping of the inflaton field begins to occur. What happens next 
is observed to depend upon the location on the k axis of the first instability band. 

1. If fc^ < for a significant fraction of the band and fc^ > for a significant fraction then 
the infiaton field goes directly into the rapid damping phase which is discussed next. 
A good example of this is the case g'^cfl = 10 which is shown in Fig. |H The average 
value of g'^{ip'^) increases exponentially until it gets large enough to significantly affect 
the infiaton field and cause the first phase of rapid damping to occur. Thus there is no 
significant amount of gradual damping that occurs before the rapid damping phase. 

2. If /c^ > for most or all of the band then as gefr decreases, the shift in the band causes 
the modes near the center to stop increasing rapidly in amplitude. In that case the 
exponential increase in {ip'^)ren ceases. At this point there are two possibilities that 
have been observed: 

(a) If the decrease in ges has resulted in fc^ > for part or all of an instability band 
in the (a, q) plane for which previously /c^ < 0, then the modes which are now in 

^ Of course, if initially fc^ > for only a small part of the first instability band, then the modes in this band 
may make a comparable or even smaller contribution to (■0^) than the modes in the second instability 
band at early times. 
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that band will begin to increase in amplitude exponentially and will do so at a 
faster rate than modes in what was previously the first instability band. 

(b) If the decrease in geff has not been large enough for a new instability band to 
appear, then slow damping will start to occur due to parametric amplification 
of those modes that are now near the center of the first instability band but 
which were not close to it before. As they become large enough to have an effect, 
Qes will gradually decrease and the slow damping which occurs is of the type 
discussed for the g'^(j)l < 1 case. Eventually Qcs will become small enough that 
an instability band that previously had A;^ < will now encompass = and 
parametric amplification will occur, again at a faster rate than that which is 
occurring in what was previously the first instability band. The modes will grow 
exponentially until rapid damping begins. This is what happens for g^^^Q = 35. 
A careful examination of the first plot in Fig. 1 shows the gradual damping phase 
followed by the rapid damping phase. 



In the rapid damping phase, as mentioned previously, KLS point out that a Minkowski 
spacetime approximation is valid and they suggest that the value of should decrease until 
it reaches about 1/4, at which point there are only narrow instability bands, so the rapid 
growth of 5'^('0^)ren ccascs, as does the fast damping of the inflaton field. They predict that 
the ratio of the amplitude of the inflaton fleld after the rapid damping to that before the 
rapid damping is 



with A the amplitude of the inflaton fleld (f) after the rapid damping, and the amplitude 
just before the rapid damping begins. We consistently find that after the rapid damping 
has finished gcfr ~ 10^^. This is significantly smaller than the value of 1/4 that they predict. 
As a result, the amount of damping which the infiaton field undergoes is much larger than 
the amount predicted by KLS. After the rapid damping, the value of qes is observed to vary 
periodically as the envelope and frequency of the oscillations of the inflaton fleld change in 
time. 



b. Rapid Damping Phase 




(4.11) 
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A careful examination of the plots showing the evolution of indicates that the rapid 
damping always seems to occur in two different phases separated by a short time. After the 
first phase the field is damped by a significant amount, but gefr is still large enough that 
significant parametric amplification can occur for modes in the lowest instability band. As 
it is occurring, but before it gets large enough to have a significant effect on g'^^il)'^)^^^, the 
inflaton field does not undergo any more noticeable damping. Once there is enough of an 
effect, the second phase of rapid damping takes place and after it is over, ges becomes small 
enough so that no more significant damping appears to occur even after a large number of 
oscillations of the infiaton field. 

It is during this period, begun by the first rapid damping of the infiaton field and finished 
by the second, that the extreme sensitivity to initial conditions that we have found seems 
to manifest. A useful way to study the sensitivity is to look at the time evolution of the 
quantities 0, g'^{il''^)ren and fk for two cases in which the initial data are almost but not 
quite identical. For the cases discussed here, 0o is different by a small amount, which then 
generates small changes in the initial values of the modes fk and thus in g'^{ip'^)rea_- If the 
relative differences in the quantities: 50/0, 5('0^)ren/('0^)ren and Sfk/fk are plotted as in 
Figs. [8] and [9l then it is seen that an exponential increase in their amplitudes takes place 
between the time when 5'^('0^)ren ~ 2, which is when the rapid damping begins, and the time 
when the average of g'^{ip'^)ren reaches its maximum value, which is when the rapid damping 
ends. 

The mechanism for this instability is not completely understood. However, from Figs. M 
and [9] it can be seen that the exponential increase in the average value of |50/0| is tied to 
that of the average value of |5('0^)/('0^)|. Further there is evidence that the average value 
of \6fk/fk\ undergoes a similar exponential increase during the same time period even for 
modes which have never undergone parametric amplification. Finally it is clear from the 
plots in Figs. HI [HI and [9] that the exponential increases occur at times when 0, ('0^)ren, and, 
at least for the mode shown, fk, are not increasing exponentially. 

To understand how this might occur, consider two cases, labeled by the subscripts 1 and 
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FIG. 8: The time evolution is shown for the case in which g = 10~^ and g'^^Q = 10 for one solution 
and 10(1 + 10~^) for the other. The top plots show the evolution of \S(f)/(l)\ and |5(V'^)ren/(^^)ren|- 
The bottom plot shows the evolution of \SRe /fc/(Re fk)\ for k = 0.5. The dashed curve corresponds 
to |Re/fc| and is not a relative difference. Note that the growth of the relative difference for each 
quantity shuts off before its average becomes of order unity. The regions of the plots where the 
average evolves as a straight line with a nonzero slope correspond to times when the average relative 
differences are growing exponentially. 

2 such that 50 = 0i — 02- Then the equations for the evolution of 50 and 5fk are 

64> = -(l + /(A')ren)50-/5(^')rcn02 , (4.12a) 

^ + 0?) Sh + 2 ihh g' 01 50 . (4.12b) 

Thus, it is clear that 5('0^) is a source for 50 and 50 is a source for Sfk- In this way, one 
can understand how the exponential increases could be correlated. 

It is not clear what is driving the exponential increase or what stops it. However, it is 
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FIG. 9: The time evolution is shown for the case in which g = 10~^ and g^cp^ = 10 for one solution 
and 10(1 + 10"^) for the other. The top plots show the evolution of \6(t)/(f>\, |(5(V'^)ren/ (^^)ren|- The 
bottom plot shows the evolution of \SRe fk / (Re fk)\ for k = 0.5. The dashed curve corresponds 
to |Re/fc| and is not a relative difference. The regions of the plots where the average evolves as 
a straight line with a nonzero slope correspond to times when the average relative differences are 
growing exponentially. 

clear that it starts when backreaction effects from the quantum field on the inflaton field 
become very important, which is just before the rapid damping occurs. It is also clear that 
it stops once the rapid damping is finished. At this point gefr <S 1. 

This appears to be a different type of sensitivity to initial conditions than that found by 
KLS. Theirs is related to sensitivity of the evolution of the phase of cf) and the value of the 
exponential parameter /j, to absolute changes in the parameter q which are of order unity, 
but which correspond to small relative changes in q for large values of q. The sensitivity 
discussed here has been observed to occur for 0.5 < g < 25. The only larger value of q that 
was investigated was q — 10^. For g = 10^ a sensitivity to initial conditions was observed 
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which becomes important at earher times than the ones we found for smaller values of q. 
Thus it seems likely that this is the sensitivity discussed by KLS. The sensitivity to initial 
conditions that we have found results in significant changes in the evolution of the inflaton 
field for changes in q which in some cases have been observed to be less than one part in 
10®, for values of q of order unity. 



c. Post-Rapid Damping Phase 



As mentioned in the Introduction, the model considered here does not include quantum 
fluctuations of the inflaton field and also does not include interactions between the created 
particles or between those particles and quantum fluctuations of the inflaton field. It would 
naturally be expected that scattering due to these interactions would become important at 



some point and that they would lead to the thermal 
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ization of the particles. Various studies 

m 
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35( 1 indicate that scattering 



which have taken such interactions into account 
effects will become important at intermediate times. In some cases this may occur before 
the rapid damping phase has ended 
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27j. Once the scattering becomes important the 



evolution should change significantly from that found using our model. 

Nevertheless it is of some interest to see what happens in our model at late times because 
of the insight it gives into the effects of particle production in the context of the semiclassical 
approximation. First, since quantum fluctuations of the inflaton field are ignored, one would 
expect that the damping of the inflaton field is due to the production of particles of the 
quantum fields. However, in the semiclassical approximation when the classical field is 
rapidly varying, different definitions of particle number can give different answers 52| . Even 
if one has a good definition, the particles may not behave like ordinary particles. This is 
the case here as can be seen in the right hand plot of Fig. [3l Recall that in Minkowski 
spacetime the total energy density is conserved. Then it is clear from the plot that after the 
rapid damping phase, most of the energy is transferred back and forth between the inflaton 
field and the quantum field. Thus, the created particles have an energy density that is very 
different from what they would have if the inflaton field was slowly varying. 

Also as mentioned previously, damping of the inflaton fleld is observed to be negligible 
for a large number of oscillations after the rapid damping phase has flnished. However, 
one would expect that since the quantum flelds are massless, particle production should 
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FIG. 10: Shown is the real part of two mode functions fk for the case g = 10 ^, g'^(j)Q = 10. For 
the left plot k = 10~^ and for the right plot k = 0.5. 



continue to occur 



531]. Examination of the behaviors of modes with various frequencies 



shows that there are some modes which do continue to grow during the post-rapid damping 
phase. One would expect that these modes would have k"^ ~ g'^A^/u^^ = 4gefr- The left 
hand plot of Fig. [10] shows that such growth does occur. However, it does not appear to 
be exponential and this may be the reason that these modes have not grown enough to 
contribute significantly to g'^{ilj'^)ren, even after a large number of oscillations of the infiaton 
field. For other modes, such as that in the plot on the right of Fig. [TOl significant growth 
ceases after the rapid damping of the infiaton field. Finally, as would be expected for modes 
with large values of A;, no significant growth in the amplitude of the oscillations was observed 
to occur. 



V. DISCUSSION AND CONCLUSIONS 

We have investigated preheating in chaotic inflation using the large expansion for a 
classical infiaton field coupled to identical massless quantum scalar fields in a Minkowski 
spacetime background. The infiaton field is minimally coupled to the gravitational field. 

When backreaction effects on the infiaton field are ignored the mode equations for the 
quantum fields are Mathieu equations. We have shown that, even though parametric am- 
plification occurs in an infinite number of bands of modes with arbitrarily large frequencies, 
at any finite time quantities such as (V'^)ren and {T^u)ren are finite due to the fact that the 
widths of the bands become negligible at very high frequencies and the maximum value of 
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the exponential parameter // in a band becomes very small at very high frequencies. 

We have investigated the backreaction effects of the quantum fields on the inflaton field 
in detail by numerically solving the full set of equations that couple the inflaton fleld to 
the quantum field. These are the first such solutions to these equations for the model we 
are considering. We have found that, in agreement with the prediction by KLS, there is a 
period of rapid damping of the infiaton field if g^^^o ~ 2 and none otherwise. In both cases 
the details of the damping can be understood through an analysis similar to that done by 
KLS; this analysis involves studying the way in which the instability bands change as the 
infiaton field damps. For g'^4>l > 2 the situation depends upon the locations of the instabihty 
bands, and therefore, there arc difi^crcnt ways in which the damping proceeds. For g'^(j)1 < 1 
there is a brief period in which the amplitude of the infiaton field is rapidly damped but 
only by a very small amount. This is followed by a long period of slow damping. Thus, we 
distinguish this from the cases in which the infiaton field is rapidly damped by a significant 
amount. 

The analysis by KLS shows that significant differences in the evolution occur for reahstic 
values of the initial conditions and parameters when the expansion of the universe is taken 
into account. However, in cases where rapid damping of the infiaton field occurs, it is possible 
to neglect this expansion during the rapid damping phase. We find that during the rapid 
damping phase there is more damping of the infiaton field than was predicted originally by 
KLS. We also find that the rapid damping occurs in two bursts separated by a period of time 
which is generally longer than the time it takes for each period of rapid damping to occur. 
Finally, we find that there is an extreme sensitivity to initial conditions which manifests in 
significant changes in the evolution of the infiaton field and the modes of the quantum fields 
during the period of rapid damping. This can occur for changes in the initial amplitude of 
the infiaton field as small as, or even smaller than, one part in 10^. Since this sensitivity 
occurs for values of g^ct)^ ^ 2, it would seem to be different than the one discussed by KLS, 
which only occurs for large values of g'^^o- 
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APPENDIX A: RENORMALIZATION AND CONSERVATION OF THE 
ENERGY-MOMENTUM TENSOR 

In this appendix the renormalization and conservation of the energy-momentum tensor 
for the system we are considering are discussed. As discussed in Section II, at leading order 
in a large N expansion, the system is equivalent to a massive classical scalar field coupled 
to a single massless quantum field. That in turn is equivalent to the system discussed in 



Ref. [47[, if the scale factor, a{t), is set equal to one, the coupling of the infiaton field to the 
scalar curvature is set to zero, the coupling constant A is set equal to 2g'^, the mass of the 
quantum field is set to zero, and the X(j)^/3\ term in the equation of motion for the field 
is dropped. Because of this, we will use several results from Ref. in what follows. 

We start with the quantum field theory in terms of bare parameters. The bare equations 
of motion for the homogeneous classical infiaton field and the quantum field modes are given 
in Eqs. (12. Sp and ( 12. 5p and an expression for the quantity {iP'^)b is given in Eq. ( 12. lip . It 
is assumed throughout that (i) the infiaton field is homogeneous and thus, depends only 
on the time t and that (ii) the quantum field t/^ is in a homogeneous and isotropic state so 
that {ip'^) depends only on time and (T^^^) when expressed in terms of the usual Cartesian 
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components, is diagonal and depends only on time. 

The total bare energy-momentum tensor of the system is given by 



47| 



rpB _ rpC,B I /rp \Q,B 



m. 



1 



2 a2 



+ 2CBV(.Ai^oij) 



(Ala) 
(Alb) 

(Ale) 



We first show that the total energy-momentum tensor is covariantly conserved. Since it is 
diagonal and depends only on time, the conservation equation in Minkowski spacetime is 



trivially satisfied for all of the components of T^, except for T^^. This component has the 



explicit form 



(Tu) 



Q,B 



ATI' 



+ 00 



dA; 



The conservation condition is 



(A2a) 
(A2b) 



B 

t J-tt 



. 



(A3) 



It is easy to show, making use of the equations of motion for (pit) and /fc(t), that 



dt {Tu) 



Q,B 



(A4a) 
(A4b) 



with {ip"^)^ given by Eq. (12. lip . This implies that the total bare energy-momentum tensor 
is covariantly conserved. 
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'. lenormalization is accomplished through the method of adiabatic regularization 



40 



41 



431]. The details are given in Ref. 47|]. The result for (V'^)b is 



(^')b - (^')ad 



27r2 
9 





2^2 



(A5a) 



(A5b) 
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For the energy-momentum tensor 



with 



{Tu)\ 



c 



i -2 m" 

+ 00 



(A6a) 



(A6b) 



47r2 



dA; A;^ 



47r2 



dkkUk + 



\fk\' + {k' + gY)\fk\ 

2 J,2" 



2A; 



1 



(Tu) 



Q 

an 



47r2 

9^ 
327r2 



,2 J,2 



2A; 



,4 J,4 



(A6c) 
(A6d) 



We note that in Minkowski spacetime the bare and the adiabatic energy density of the 
quantum field do not depend on the value of the coupling ^ of the quantum field to the 
scalar curvature. Thus, the renormalized value of the energy density does not depend on ^. 

Using Eqs. (lA6p . (lASp . and the equations of motion for and /fc, it is easy to derive the 
following identities: 



dtT^t 



They can be used to show that 

dt {Tu)ren = dt {Tttf + Ot {Ttt)i " dt {Tttf,^ 

= [-(^')ren+(V^^)B-(V^^)ad] 



(A7a) 
(A7b) 
(A7c) 



(A8) 
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